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Abstract. We obtain the hydrodynamic limit of a simple exclusion process in an 
inhomogeneous environment of divergence form. Our main assumption is a suitable 
version of T-convergence for the environment. In this way we obtain an unified 
approach to recent works on the field. 



1. Introduction 

Since the seminal paper Guo et al. (1988), the theory of hydrodynamic limit of 
interacting particle systems has evolved into a powerful tool in the study of non- 
equilibrium properties of statistical systems of many components (see the book 
Kipnis and Landim (1999) for a comprehensive exposition). Recently, and due to 
the infuence of physical and mathematical works about random walks in random 
environment, an increasing attention has been posed into particle systems evolving 
in random environments. Despite the early works Fritz (1989), Quastel (2006), 
Koukkous (1999), we mention Faggionato and Martinelli (2003), Quastel (2006) 
Nagy (2002), Jara and Landim (2006), Faggionato (2007), Gongalves and Jara 
(2008a), Faggionato et al. (2009), Franco and Landim (2008), Faggionato (2008), 
Gongalves and Jara (2008b). In Gongalves and Jara (2008a), Jara and Landim 
(2006) the corrected empirical density was introduced, which is nothing but a mi- 
croscopic version of the compensated compactness lemma of Tartar Tartar (1979). 
Roughly speaking, when the inhomogeneous environment (random or not) has a 
divergence form and has a T- limit, space homogenization of the environment and 
time homogenization of the interaction decouples, and the standard tools from the 
theory of hydrodynamic limit can be used to obtain the asymptotic behavior of the 
density of particles in a family of models, including the exclusion process and the 
zero-range process. 

In this review, we give an unified approach to this problem, recovering previous 
results in Nagy (2002), Jara and Landim (2006), Faggionato (2007), Faggionato 
et al. (2009), Faggionato (2008), in a simple way. In order to concentrate our efforts 
in the influence of the inhomogeneous environment on the asymptotics of the density 
of particles, we consider the simplest model of interacting particle systems, which 
is the symmetric exclusion process 77" in an unoricntcd graph. In this process, 
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particles perform symmetric random walks on a graph {X n } n with some rates 
uj n = {ui™ y ;x,y £ X n }, conditioned to have at most one particle per site. Wc 
think of {X n } n as a sequence of graphs embedding in some metric space X, and 
we are interested in the evolution of the measure 7r™(dx) in X, obtained by giving 
a mass a^ 1 to each particle. 

This article is organized as follows. In Section 2 we give precise definitions 
of the exclusion process, the inhomogeneous environment and we state our main 
result. We also define what we mean by an approximation {X n } n of X and by F- 
convegence of the environment. In Section 3 we introduce the corrected empirical 
density and we prove our main theorem. In Section 4 we introduce the concept 
of energy solutions of the hydrodynamic equation, we prove uniqueness of such 
solutions and we obtain a substantial improvement of the main Theorem. The 
material of this Section is new and it gives a better understanding of the relation 
between T-convergence of the environment and hydrodynamic limit of the particle 
system. In Section 5 we discuss how to rcobtain previous results in the literature 
relying in our main Theorem. 

2. Definitions and results 

In this section we define the exclusion process in inhomogeneous environment 
and we recall some notions of T-convergence that will be necessary in order to 
obtain the hydrodynamic limit of this process. 

2.1. Partitions of the unity and approximating sequences. In this section we fix 
some notation and we define some objects which will be useful in the sequel. Let 
(X, B) be a Polish space. We assume that X is cr-compact. We say that a sequence 
of functions {U^ i £ 1} is a partition of the unity if: 

i) : for any i £ 7, Ui : X — » [0, 1] is a continuous function, 

ii) : for any x £ X, Y.iei u i( x ) = ^ 

iii) : for any x £ X, the set {i £ I;Ui(x) > 0} is finite. 

We say that the partition of the unity {Ui ; i £ 1} is regular if supp Ui is compact 
for any i £ 7, and additionally Ui(X) = [0,1]. We denote by A4 + (X) the set 
of Radon, positive measures in X. The symbol {x n } n will denote a sequence of 
elements x n in some space, indexed by the set N of positive integers. 

Let {Ui}i be a regular partition of the unity. We say that a sequence {xi\ i £ 1} 
in X is a representative of {Ui}i if Ui(xi) = 1 for any i £ I. Notice that we have 
Xi 7^ Xj for i ^ j. 

Let {UI 1 ; i £ I n }n be a sequence of regular partitions of the unity. We say that 
a measure (i £ A4+(X) is the scaling limit of the sequence {U™} n if there exists 
a sequence {a„}„ of positive numbers such that for any sequence {x™;i £ I n } of 
representatives of {U™} n we have 



with respect to the vague topology, where S x is the Dirac mass at x £ X. We call 
{a n }n the scaling sequence. 



lim — 

rwoo a n 
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From now on, we fix a sequence of regular partitions of the unity with 

scaling limit /j,, scaling sequence {a„}„ and we assume that fi(A) > for any non- 
empty, open set AC X. Fix a sequence {xf;i g I n } of representatives of 
Define X n = {xf;i g /„}. Since {U™} is a partition of the unity, the induced 
topology in X n coincides with the discrete topology. For i = i", we will denote 
W^U^. Define 

Hn(dx) = — ^2 S x (dx). 

x£X n 

By definition, fi n — > \i in the vague topology. We denote by £ 2 (yU„) the Hilbert 
space of functions / : X n — > R such that X^ex /( x ) 2 < +0°, equipped with the 
inner product 

(/,ff>n = 7- £ f{x)g{x). 

n xex n 

We define £ 2 (p) , £ : (A„) and C 1 (/z) in the analogous way and we denote (/, g) = 
j fgdfi. We denote by C C (X) the set of continuous functions / : X — > R with 
compact support. In the same spirit, we denote by C C (A„) the set of functions 
/ : X n — ► R with finite support. We define the projection £„ : C C (X) — » C C (A„) by 
taking 



(S n G){x) = a n f GUSdfi. 



This operator, under suitable conditions, can be extended to a bounded oper- 
ator from C 2 (X) to C 2 (X n ). Notice that J S n Gdfj, n = j Gdfj,. Therefore S n is 
continuous from to £ 1 (X„). 

2.2. T -convergence. Define R = [—00, +00]. Let (Y,F) be a topological space, and 
let F n , F : Y -> R. We say that F„ is T-convergent to F if: 

i) : For any sequence {y n }n in Y converging to y g F, 

F(y) <lim inf F„ (?/„). 

ii) : For any y g y there exists a sequence {j/„}„ converging to y such that 

limsupF„(j/„) < F(y). 

n — >oo 

An important property of T-convergence is that it implies convergence of mini- 
mizers in the following sense: 

Proposition 2.1. Let F n , F : Y — > R &e s«c/i i/i<z£ F„ is Y -convergent to F . Assume 
that there exists a relatively compact set K C y sucA i/iai /or any n 7 



TTien, 



inf F„(i/)= inf F„(y). 



lim inf F n (y) = min F(y). 

n—>ocy£K " y£Y 



Moreover, if {y n } n is a sequence in K such that lim n (F n (y n ) — mipc F n ) — 0, then 
any limit point y of {y n } n satisfies F(y) = miny F. 

A useful property that follows easily from the definition is the stability of T- 
convergence under continuous perturbations: 
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Proposition 2.2. Let F n ,F : Y — > R be such that F n is T-convergent to F. Let 
G n : Y — > R be such that G n converges uniformly to a continuous limit G. Then, 
F n + G n is r -convergent to F + G. 

2.3. The exclusion process in inhomogeneous environment. In this section we de- 
fine the exclusion process in inhomogeneous environment as a system of particles 
evolving in the set X n . Let ui n = {uj™ y ;x,y G X n } be a sequence of non-negative 
numbers such that w™ x = and = w™ x for any x,y G X„. We call u n 
the environment. We define the exclusion process 77™ with environment u™ as a 
continuous-time Markov chain of state space 0„ = {0, l} Xn and generated by the 
operator 

L n f( V )= £ < y [f(v x ' y )-m], 

x,y£X n 

where 77 is a generic element of fl n , f : fl n — > R is a function which depends on 77(2;) 
for a finite number of elements x G X n (that is, / is a local function) and ry 2 ' 2 ' G £1„ 
is defined by 

iiz = x 
n(x), if z = y 
T](z), if zj^x,y. 

In order to have a well-defined Markovian evolution for any initial distribution 
t/q , we assume that sup^, ^2 yGXn u x, v < +°°- We interpret X n as a set of sites and 
rjt(x) as the number of particles at site x G X n at time t. Since r/™(a;) e {0,1}, 
there is at most one particle per site at any given time: this is the so-called exclusion 
rule. Notice that the dynamics is conservative in the sense that no particles are 
annihilated or destroyed. 

Our interest is to study the collective behavior of particles for the sequence of 
processes {n n } n - In order to do this, we introduce the empirical density of particles 
as the measure- valued process 7r™ defined by 

7r t "(G) = -L ]T r£{x)S n G{x) 

a n v 

for any G G C C (X). Using Riesz's theorem, it is not difficult to check that 7r™ 
is effectively a positive Radon measure in X. Observe that when vft{x) = 1 for 
any x G X n , then 77™ (x) = 1 for any x G X n and any t > 0. In this situation, 
the empirical process 7r™ is identically equal to the measure /1. Notice that the 
random variable 7r™ defined in this way corresponds to a process defined in the 
space T>([0, 00), A4+(X)) of cadlag paths with values in M + {X). For functions 
G : X n -> R, we define tt?(G) = a" 1 rj?(x)G(x). 

2.4. Y -convergence of the environment. In this section we will make a set of as- 
sumptions on the environment {w™}„ which will allows us to obtain an asymptotic 
result for the sequence {7r™}„. We start with two assumptions about the sequence 
of partitions of the unity Our first assumption corresponds to a sort of 
ellipticity condition on the partitions of the unity {U"} n - 

(HI): There exists O < +00 such that 

sup a n / U^dfi < 6 for any n > 0. 
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Under this condition, the projection S n satisfies ||S n G||oo < #||G||oo, and by 
interpolation S n can be extended to a continuous operator from £ 2 (/x) to £ 2 (X n ). 
Our second condition states that S n is close to an isometry when n — > oo: 

(H2): For any F g £ 2 (m), we have 

lim (S n F,S n F) n = (F,F). 

Now we are ready to discuss on which sense we will say that the environment u n 
converges. For a given function F : X n — ► R of finite support, we define £ n F by 

C n F{x)= u>Z v (F{y)-F(x)). 
yex n 

It turns out that C n can be extended to a non-positive operator in C 2 {X n ). In 
fact, for any function F of finite support, the Dirichlet form 

(F,-C n F) n = ±- ]T ul y (F(y)-F{x)) 2 

is clearly non-negative. For a function G G £ 2 (^), define £ n (G) = (S n G, —C n S n G). 
Notice that £ n : C 2 (ji) — > R is a quadratic form. Now we are ready to state our 
first hypothesis about the environment: 

(H3): There exists a non- negative, symmetric operator C : D(C) C £ 2 (yu) — > 
£ 2 (^t) such that £„ is T-convergent to £, where £{G) = — J GLGdfi. 

Our second hypothesis about the environment u n concerns to its T-limit C: 

(H4): There exists a dense set JC C C C (X) such that K. is a kernel for the 
operator £, and for any G e /C, £G is continuous and J |£G|rf^i < +00. 



2.5. Hydrodynamic limit of rff . In this section we explain what we understand as 
the hydrodynamic limit of rj". We say that a sequence {is n } n of distributions in f2„ 
is associated to a function u : X — > R if for any function G € C C (X ) and any e > 
we have 

lim f„{|— rj(x)G(x) — / G(x)u(x)fi(dx) > e\ = 0. 

Notice that we necessarily have < u(x) < 1 for any x e X, since rj(x) g {0, 1}. 
Fix an initial profile m : -X" — ► [0, 1] and take a sequence of distributions {v n } 
associated to uq. Let rff be the exclusion process with initial distribution v n . 
We denote by P„ the law of 77™ in T>([0, 00), Q n ) and by E„ the expectation with 
respect to P„. The fact that {v n } n is associated to uo can be interpreted as a law 
of large numbers for the empirical measure tTq : ttq (dx) converges in probability to 
the deterministic measure uo(x) fj,(dx) . We say that the hydrodynamic limit of jj" 
is given by the equation d t u = Cu if for any t > 0, the empirical measure w"(dx) 
converges in probability to the measure u(t,x)fj,(dx), where u(t,x) is the solution 
of the equation d t u = Cu with initial condition u . Before stating our main result 
in a more precise way, we need some definitions. 

For F,G g D{C), define the bilinear form £(F,G) = - J FCGd/j,. Notice that 
£{F, G) is still well defined if only G g D(C). We say that a function u : [0, T] xX -» 
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[0, 1] is a weak solution of (2.1) with initial condition uq if J J u\d[idt < +oo and 
for any differentiable path G : [0, T] — > )C such that Gt = we have 



Theorem 2.3. Let {v n }„ be associated to uq and consider the exclusion process 
rft with initial distribution v n . Assume that J ^(dx) is uniformly finite: 



Then, the sequence of processes {Tr n (dx)} n is tight and the limit points are con- 
centrated on measures of the form u{t,x)ji{dx), where u(t,x) is a weak solution of 
the hydrodynamic equation 



If such solution is unique, the process ir n (dx) converges in probability with re- 
spect to the Skorohod topology ofT>([0, oo), M.+(X)) to the deterministic trajectory 
u(t, x)fi(dx). 

Usually in the literature, hydrodynamic limits are obtained in finite volume, 
since the pass from finite to infinite volume is non-trivial. Assumption (H5) is in 
this spirit: it is automatically satisfied when the cardinality of X n is of the order 
of a n (on which case n{X) < +oo), and it is very restrictive when X n is infinite. 
For simplicity, we restrict ourselves to the case on which (H5) is satisfied. 

3. Hydrodynamic limit of n": proofs 

In this section we obtain the hydrodynamic limit of the process 77™ . The strategy 
of proof of this result is the usual one for convergence of stochastic processes. First 
we prove tightness of the sequence of processes {7r™}„. Then we prove that any limit 
point of this sequence is concentrated on solutions of the hydrodynamic equation. 
Finally, a uniqueness result for such solutions allows us to conclude the proof. 
However, the strategy outlined above will not be carried out for {ir n } n directly, but 
for another process 7r™, which we call the corrected empirical process. 

3.1. The corrected empirical measure. In this section we define the so-called cor- 
rected empirical measure, relying on the T-convergence of the environment. First 
we need to extract some information about convergence of the operators C n to C 
from the T-convergence of the associated Dirichlet forms. 

Take a general Hilbert space H and let A be a non-negative, symmetric operator 
defined in H. By Lax-Milgram theorem, we know that for any A > and any g eH, 
the equation (A + A)f = g has a unique solution in TL. Moreover, the solution / is 
the minimizcr of the functional / ^ (f,Af) + A||/|| 2 — 2(f,g). Fix A > 0. For a 
given function G G £ 2 (m), define the functional 




(H5): 




dtu — Cu, 
u(0, •) = m o (0- 



(2.1) 




(F)=£ n (F) + X(S n F,S n F) r . 
£ G (F)=£(F) + X{F,F) 



n 



2(S n F, S, 



•n 



G) 



2(F,G). 
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By Proposition 2.2, £ G is T-convergent to £ G . In particular, a sequence of 
minimizers F n of £ G converge to the minimizer F of £ G . Notice that F n is not 
uniquely defined in general, although S n F n it is. By the discussion above, (A — 
£ n )S n F n = S n G and (A — C)F = G. Since the operator norm of S n is bounded by 
6, we conclude that the £ 2 (A„)-norm of S n F n — S n F converges to as n — > oo. 
By (H2), we conclude that £ n (F n ) converges to £(F). 

Now we are ready to define the corrected empirical measure 7r™ . Take a function 
Ge/C and define H = (A — C)G. Define G n as a minimizer of £^ ■ Notice that in 
this way S n G n is uniquely defined. Then we define 

K(G) = — £ V?(x)S n G n (x). 

x£X n 

In order to prove that 7r"(G) is well defined, we need to prove that J2x S n G n (x) 
is finite. Remember that (A — C n )S n G n = S n H. Consider the continuous-time 
random walk with jump rates y . Remember that the condition sup^ ~^2 y y en- 
sures that this random walk is well defined. Let y) be its transition probability 
function. An explicit formula for S n G n in terms of p^(x,y) is 

poo 

S n G n {x) = / e~ At V p^{x,y)S n H{y)dt. 
Jo yex n 

Since ^2 x Pt(x,y) = 1 for any y e X t , we conclude that 



- V S n G n (x) = \ [HdfjL 
ln xex n A J 



and in particular S n G n is summable. We conclude that 7r"(C?) is well defined. 
Notice that it is not clear at all if 7r" is well defined as a measure in X. 

3.2. Tightness of {n n } n and proof of Theorem 2.3. In this section we prove tight- 
ness of {TT n } n and we prove Theorem 2.3. As we will see, we rely on the corrected 
empirical measure, which turns out to be the right object to be studied. By (H5), 
we have 

lim F n ( sup k n (G) - 7r t "(G)| > e) = 0. 

Notice that (H5) can be substituted by the following condition, which can be 
sometimes proved directly. 
(H5')= ForanyGe/C, 



lim — V \S n G n (x) - S n G(x)\ =0. 



n^oc On 



In particular, {ir n (G)} n is tight if and only if {7r™(G)}„ is tight. The usual way 
of proving tightness of {n n (G)} n is to use a proper martingale decomposition. A 
simple computation based on Dynkin's formula shows that 

M n t {G) = 7r t "(G) - ^ l (G) - /' n n s {C n S n G n )ds (3.1) 

Jo 

is a martingale. The quadratic variation otM"(G) is given by 

(MUG))=f t ^ E (r£{v)- V ?{x)) 2 wZ y (S n G n {v)-S n G n {x)) 2 ds. 
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In particular, (A4f(G)) < ta~ 1 £ n (G n ). At this point, the convenience of intro- 
ducing the corrected empirical process becomes evident. By definition, C n S n G n — 
S n CG + \(S n G n — S n G). Since H = (A — C)G, the function G is the minimizcr 
of £ H . Therefore, G n converges to G in C 2 (X). By (H2), the £ 2 (X„)-norm of 
S n G n — S n G goes to and £ n (G n ) converges to £(G). 

We conclude that A4™(G) converges to as n — > oo, and in particular the 
sequence {M n (G)} n is tight. In the other hand, the integral term in (3.1) is equal 
to f*w?(£G)ds. 

Notice that 7r"(£G) < / \£G\d[i for any t > 0, from where we conclude that the 
integral term is of bounded variation, uniformly in n. Tightness follows at once. 
Since {ttq (G)}„ is tight by comparison with {ttq (G)}„, we conclude that {7r"(G)}„ 
is tight, which proves the first part of Theorem 2.3. As a by-product, we have 
obtained tightness for {n n } n as well, and the convergence result 

lim L»(G) - tt "(G) - (\ n s {LG)ds\ = 

for any G e JC. Notice that we have exchanged tt"(G) by 7r™(G). Let ir. be a limit 
point of {n n } n . Then, ir. satisfies the identity 

7r t (G)-7r (G)- / 7r s (£G)ds = 
Jo 

for any function G € JC. By hypothesis, TTo(dx) = uo(x)(i(dx). Repeating the 
arguments for a function G t (x) = Gq(x) + tG\{x) with Go, Gi G /C, we can prove 
that 

7T t (G t )-n (G )- [ ir s ((d t +£)G s )ds = 
Jo 

for any piecewise-linear trajectory G. : [0, T] — ► JC. The same identity holds by 
approximation for any smooth path G. : [0,T] — > C C (X), which proves that the 
process it. is concentrated on weak solutions of the hydrodynamic equation. When 
such solutions are unique, the process ir is just a ^-distribution concentrated on 
the path u(t,x)fi(dx). Since compactness plus uniqueness of limit points imply 
convergence, Theorem 2.3 is proved. 

4. Energy solutions and energy estimate 

In this section we define what we mean by energy solutions of Equation (2.1), we 
prove that any limit point of the empirical measure {ir n } is concentrated on energy 
solutions of (2.1) and we give a simple criterion for uniqueness of such solutions. 

4.1. Energy solutions. Let £ : H — > M be a quadratic form defined over a Hilbcrt 
space H of inner product (•, •). We say that £ is closable if for any sequence {/„}„ 
converging in H to some limit / such that £( f n ~ fm) goes to as n, m — > oo, we 
have / = 0. Let £ : H — > M be closable. We define Hi = TL\{£) as the closure of 
the set {/ e H;£(f) < +00} under the norm ||/||i = (£(f) + (/, f)) 1 ' 2 . 

We say that a dense set K C H is a kernel of £ if Wi is equal to the closure of 
K under the norm || • ||i. We say that a symmetric operator C : D(C) C H — > i7 
generates 5 if £(/) = (/, — £/) for / e D(£) and D(£) is a kernel of £. 
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Fix T > 0. For a function u : [0, T] — > ii we define the norm 

I KIM*) 

and we define Wi,t as the Hilbert space generated by this norm. Given a closable 
form £ generated by the operator £, we say that a trajectory u : [0, T] — > H 
is an energy solution of (2.1) if u € Wi,t and for any diffcrcntiable trajectory 
G : [0, T] -> Hi with G(T) = we have 

(Go, « > + jf {(0 t G t , u t ) - £{G t ,u t ))dt = 0. 

In other words, an energy solution of (2.1) is basically a weak solution belonging 
to Wi,t- In fact, by taking suitable approximations of G, it is enough to prove this 
identity for trajectories G such that G t G K for any t G [0,T], where if is any 
kernel of £ contained in D(C). Notice that the norm in Ti\ t T is stronger than the 
norm u^dt, and therefore a weak solution is effectively weaker than an energy 
solution of (2.1). 

4.2. The energy estimate. In this section we prove that the limit points of the 
empirical measure are concentrated on energy solutions of (2.1). For simplicity, we 
work on finite volume. From now on we assume that X is compact. Therefore, 
there exists a constant k such that the cadinality of X n is bounded by na n . We 
have the following estimate. 

Theorem 4.1. Fix T > 0. Let {H l : X n x X n x [0, T] -> M; i = 1, . . . , 1} be a finite 
sequence of functions. There exists a constant C = C{T) such that 



i=l,...,lJ x . yeXn 

E <A)Vw}*l<c+^. (4.1) 

Proof. Before starting the proof of this theorem, we need some definitions. Fix 
p > 0. Denote by v p the product measure in Q n defined by 

uP(r)( Xl ) = l,...,r](x k ) = 1) =p k . 

It is not difficult to check that the measure v p is left invariant under the evolution 
of rjt- For two given probability measures P±, P2, we define the entropy H(P\\P2) 
of Pi with respect to P2 as 



H(P!\P 2 ) 



+00, if Pi is not absolutely continuous with respect to P 2 

/ log ^ rfPi otherwise . 



For rj G f2 n , denote by 5 V the Dirac measure at n. It is not difficult to see 
that H{8 v \v p ) < C(p)a n for any r\ G O n , where G(p) is a constant that can be 
chosen independently from n. Let us denote by P p the distribution in D([0, T],Sl n ) 
of the process 77™ with initial distribution By the convexity of the entropy, 
H(¥ n \FP) < C(p,T)a n for a constant C(p,T) not depending on n. The following 
arguments are standard and can be found in full rigor in Kipnis and Landim (1999). 
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Let us denote by F l (s) the function (depending on H l (s) and 77™) under the time 
integral in (4.1). By the entropy estimate, 



E„ 



sup / 

i=l,...,l JO 



F l {t)dt 



< 



H(¥ n \FP) 



+ 



1 r r T 

— logE p exp{ sup a n / F l {t)dt) 

a n L i=l,...,J Jo 



In order to take the supremum out of the expectation, we use the inequalities 
expjsupi&i} < J2i ex P{bi} an d l°g{Z)» &*} < log / + supjlog&j, valid for any real 
numbers {bi, i = 1, . . . , I}. In this way we obtain the bound 



sup I 1 F l (t)dt] < C{p, T) + ^ 



+ sup — logE p 

i=l,...,l a n 



exp{a„ / F l (t)dt} 
Jo 



(4.2) 



Therefore, it is left to prove that the last supremum is not positive. It is enough 
to prove that the expectation E p [exp { J Q T F l (t)dt}] is less or equal than 1 for any 
function F l . From now on we drop the index i. By Feynman-Kac's formula plus 
the variational formula for the largest eigenvalue of the operator F(t) + L n , we have 

— logE"[cxp{a„ f F(t)dt}]< [ sup {{F(t),f) p - {f, -L n f) p }, 
a n 1 Jo 1 Jo f 

where we have denoted by (•, -) p the inner product in C 2 {v p ) and the supremum is 
over functions / € C 2 {v p ). A simple computation using the invariance of v p shows 
that 

(f,-L n f) p = £ < y l[f{rf»)-m}\{dri). 

x,yeX n J 

Recall the expression for Fit) in terms of H . We will estimate each term of the 
form 2a- 1 (H x . v iriiy) - T](x)), f 2 ) p separatedly: 

x 

2 f {FIx,y) 2 Px,y 



-(H x , y (v(y) V(x)), f) P = —H x , v { V (x),f{rf'V) 2 fiv) 2 ) P 

Ctr), Q*r). 



< 



{ 



(v(x),(f(v x ' y )+mf) P 



+ 



Choosing (f£ = y and putting this estimate back into (4.2), we obtain the 
desired estimate. □ □ 

Take G % e JC and take H* = S n G l n iy) — S n G % n {x), with G l n defined as in Section 
3.1. Recall the identity £„5„Gj l = S n CG i + \{S n G l n - S n G' 1 ). The energy estimate 
(4.1) gives 



E rl 



sup / (27r?(rG i ) - SniG^dt <C(p,T) + C 1 il,n), 

i=l,...,lJ0 J 



where Ci(l,n) is a constant that goes to when I is fixed and n — > 00. Take a 
limit point of the sequence {tt"}„. We have already seen that 7r"(£G I ) converges 
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to Tr t (£G). Therefore, the process ir. satisfies 

sup / (2n s (£G l ) -£{G l ))dt\ < C(p,T). 
=i,...,iJo i 



Similar arguments prove that for piecewise linear trajectories {G\; i — 1, . . . , 1} 
in JC, we have 



sup / (2n s (£G i (t)) - £{G\t)))dt\ < C(p,T). 



E 

i=l, 



Since I is arbitrary and piecewise linear trajectories with values in IC are dense 
in Hi t T, we conclude that £7[| |7r. \\f T ] < +oo, from where we conclude that 1 1 7r. 1 1 1 ,r 
is finite a.s. We establish this result as a theorem. 

Theorem 4.2. Let 77™ an exclusion process as in Theorem 2.3. If one of the 
following conditions is satisfied, 

i) X is compact, 

ii) Assumption (H5') holds and the entropy density is finite: 

ff(Pn|P") , 

sup < +OO, 

n a n 

then any limit point of the sequence {ir n (dx)} n is concentrated on energy solutions 
of the hydrodynamic equation (2.1). In particular, since such energy solutions are 
unique, the sequence {-K n (dx)} n is convergent. 

4.3. Uniqueness of energy solutions. In this section we prove uniqueness of energy 
solutions for (2.1). Since the equation is linear, it is enough to prove uniqueness 
for the case uq = 0. Let u t be a solution of (2.1) with u = 0. Then, 



L 



T 

{(d t G t ,u t )-£(G t ,u t )}dt = 



for any differentiable trajectory in TL\^t with Gt = 0. Take Gt — — u s ds. Then 

dtGt — u t and the first term above is equal to (ut,u t )dt. An approximation 
procedure and Fubini's theorem shows that the second term above is equal to 

\s([u t dt). 

Both terms are non-negative, so we conclude that J Q T ' (u t , u t )dt = and u t = 0. 
5. Applications 

In this section we give some examples of systems on which Theorems 2.3 and 4.2 
apply. In the literature, the sequence uj n is often referred as the set of conductances 
of the model. Unless stated explicitely, in these examples, X will be equal to M. d 
or the torus T d = R d /Z d . The set X n will be equal to n~ 1 Z d and we construct the 
partitions {U™} in the canonical way, taking U™ as a continuous, piecewise linear 
function with U™{x) = 1 and U"(y) = for y e X n , y 7^ x. 
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5.1. Homogenization of ergodic, elliptic environments. Let P) be a proba- 
bility space. Let {t x ;x £ Z d } be a family of .F-mesurable maps t x : Q — > f2 such 
that 

i) : P^A) = P(A) for anyAef,i€ Z d . 

ii) : t x t x i = T x+X > for any x,x' £ 7L d . 

iii) : If t x A = A for any x £ Z d , then P(A) = or 1. 

In this case, we say that the family {T x } x< z%d is ergodic and invariant under P. 
Let a = (oi,.. . , ad) : £1 — > K d be an JF- measurable function. Assume that there 
exists eo > such that 

e < o»(w) < e^ 1 for all w g O and i = 1, . . . , d. 

We say in this situation that the environment satisfies the ellipticity condition. 
Fix lo £ fl. Define u> n by w" , , = uj n , , = n 2 ai(T n XLo), oj™ „ — if |y — x| ^ 

>> x,x+ei/n x+ei/n,x f V /' z.y If I / 

l/n. Here {e^ji is the canonical basis of Z d . In this case, a n — n d and \i is the 
Lebesgue measure in W d . In Papanicolaou and Varadhan (1982), it is proved that 
there is a positive definite matrix A such that the quadratic form £ n associated to 
oj n is F-convergent to £(f) = fVf ■ AVfdx, P — a.s. In particular, Theorem 2.3 
applies with Cf = div(AV/). This result was first obtained in Goncalves and Jara 
(2008a). 

5.2. The percolation cluster. Let e = {e x ;x £ Z d ,i = l,--- ,d} be a sequence of 
i.i.d. random variables, with P(e x = 1) = 1 — P(e x = 0) = p for some p = (0, 1). 
Define for x,y £ X n , u™ x+eJn = ^ +eJn x = n 2 e^, < y = if \y - x\ ± l/n. Fix 
a realization of e. We say that two points x, y £ X n are connected if there is a finite 
sequence {x = x, . . . ,xi = y} C X„ such that — orj| = l/n and w" j = 1 for 
any i. Denote by Co the set of points connected to the origin. It is well known that 
there exists p c £ (0, 1) such that 9{p) = P(Cq is infinite ) is for p < p c and positive 
for p > p c . Fix p > p c . Define a n — n d and [io(dx) = 9(p)dx. In ?, it is proved 
that there exists a constant D such that, P — a.s in the set {Co is infinite }, the 
quadratic form £ n associated to the environment uj n restricted to Co is T-convergent 
to £{f) = 9{p)D f (V f) 2 dx. Theorem 2.3 applies with C = DA, assuming that the 
initial measures v n put mass zero in configurations with particles outside Co- This 
result was first obtained in Faggionato (2007), relying on a duality representation 
of the simple exclusion process. 

5.3. One- dimensional, inhomogeneous environments. In dimension d = 1, the T- 
convergence of £ n can be studied explicitely. For nearest-neighbors environments 
(w™ y = if | a; — y\ = 1), T-convergence of £ n is equivalent to convergence in 
distribution of the measures 

W n (dx) = 1 ^{^^ x+1 )- l 5 x/n {dx). 

Let W(dx) be the limit. We assume that W(dx) gives positive mass to any open 
set. For simplicity, suppose that JF"({0}) = 0. Otherwise, we simply change the 
origin to another point with mass zero. For two functions /,j:t^Ewc say that 
g = df/dW if 

f(x) = f(0)+ f g{y)W{dy). 
Jo 
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Then £ n is r-convergent to the quadratic form defined by £ (/) = J (df /dW) 2 dW . 
In this case, C — d/dxd/dW. A technical difficulty appears if W(dx) has atoms. 
In that case, there is no kernel K, for C contained in C C (R). To overcome this 
point, we define for x < y, dw{x,y) — dw{y,x) = W((x,y]). The function dw is 
a metric in R, and in general R is not complete under this metric: an increasing 
sequence x n converging to x is always a Cauchy sequence with respect to dw, 
but dw(x n ,x) > W({x}), which is non-zero if x is an atom of W. Define Ru/ = 
R U {x— ; VF({x}) > 0}. It is easy to see that Rw/ is a complete, separable space 
under the natural extension of dyy, and that continuous functions in W.\y are in 
bijection with cadlag functions in R with discontinuity points contained on the set 
of atoms of W(dx). It is not difficult to see that the set of IF-diffcrentiable functions 
in C c (K.w) is a kernel for C and that Theorems 2.3 and 4.2 apply to this setting. 
In Faggionato et al. (2009), the remarkable case on which W(dx) is a random, 
self-similar measure (an a-stable subordinator) was studied in great detail. 

5.4. Finitely ramified fractals. Let us consider the following sequence of graphs in 
R 2 . Define a = (0,0), a 1 = (1/2, V3/2), and a 2 = (1,0) and define ifi : R 2 -» R 2 
by taking ipi(x) = (x + Oj)/2. Define X a = {00,0,02} and X n+ i = Uj<£,(X n ) for 
n > 0. For x,y e X we define = 1, we put u® y — if {x,y} C X and 
inductively we define 

, ,n+l _ rc\^ , ,n 

The set X„ is a discrete approximation of the Sierpinski gasket X defined as the 
unique compact, non-empty set X such that X = Linfi(X). Here we are just saying 
that w™ = 5™ if x, y are neighbors in the canonical sense. In this case a n = 3™ 
and /j, is the Hausdorff measure in X. It has been proved Kigami (2001) that the 
quadratic forms £ n converge to a certain Dirichlet form £ which is used to define an 
abstract Laplacian in X . In particular, Theorems 2.3 and 4.2 apply to this model. 
This result was obtained in Jara (2009) in the context of a zero-range process. The 
same result can be proved for general finitely ramified fractals, in the framework of 
Kigami (2001). 
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